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g : Abstract 

' We report on transcritical bifurcations of periodic orbits in non-integrable two-dimensional Hamil- 

^ ■ tonian systems. We discuss their existence criteria and some of their properties using a recent math- 

\ ematical description of transcritical bifurcations in families of symplectic maps. We then present 

numerical examples of transcritical bifurcations in a class of generalized Henon-Heiles Hamiltonians 
and illustrate their stabilities and unfoldings under various perturbations of the Hamiltonians. We 
demonstrate that for Hamiltonians containing straight-line librating orbits, the transcritical bifur- 
cation of these orbits is the typical case which occurs also in the absence of any discrete symmetries, 
while their isochronous pitchfork bifurcation is an exception. We determine the normal forms of 
both types of bifurcations and derive the uniform approximation required to include transcritically 
bifurcating orbits in the semiclassical trace formula for the density of states of the quantum Hamil- 
tonian. We compute the coarse-grained density of states in a specific example both semiclassically 
and quantum mechanically and find excellent agreement of the results. 

> 

m : 1 Introduction 
in 

. The transcritical bifurcation (TCB), in which a pair of stable and unstable fixed points of a map 

exchange their stabilities, is a well-known phenomenon in one-dimensional non-Hamiltonian systems. 
O ' A simple example occurs in the quadratic logistic map (see, e.g., jl]): 
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Xn+l = r Xn{l - Xn) , (1) 



\ where are arbitrary real numbers and r is the control parameter. This map has - amongst 

^ ' others - two fixed points x\ = Q and = 1 — 1/r which exchange their stabilities at the critical value 

r = 1. For r < 1, is stable and is unstable, whereas the inverse is true for r > 1. [Note that 
in many textbooks discussing the quadratic map, this bifurcation is not mentioned as the values of 
the variable x are usually confined to be non- negative, while < for r < 1.] The TCB occurs in 
many maps used to describe growth or population phenomena (see [2j for a recent example). TCBs 
have also been reported to occur in various time-dependent model systems O HI El El El El |^ and 
shown, e.g., to be involved in synchronization mechanisms |5|l6j. In [8l|9], TCBs have been found 
to play a crucial role in transitions between low- and high-confinement states in confined plasmas, 
and their unfoldings have been analyzed. 

To our knowledge, TCBs have not been discussed so far for periodic orbits in autonomous 
Hamiltonian systems. In this paper we report on the occurrence of such bifurcations in a class of 
two-dimensional non-integrable Hamiltonian systems. Since the TCB does not belong to the generic 
bifurcations in two-dimensional symplectic maps [10], we consider it useful to investigate also the 
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mathematical conditions under which it can exist, its stabihty under perturbations of the Hamilto- 
nian, and its unfoldings when it is destroyed by a perturbation. For this, we rely on mathematical 
studies by Janich [111 [T2] who introduced a class of "crossing bifurcations", to which the TCB 
belongs, and derived several theorems and useful formulae for crossing bifurcations of straight-line 
librational orbits. Finally, in view of the important role which Gutzwiller's semiclassical trace for- 
mula ^13j plays for investigations of "quantum chaos" (see, e.g., [HI IS]), we study the inclusion of 
transcritically bifurcating orbits in the trace formula by an appropriate uniform approximation. 

Generic bifurcations of fixed points in two-dimensional symplectic maps have been classified by 
Meyer [TU] in terms of the number m = 1, 2, . . . that corresponds to a period m-tupling occurring 
at the bifurcation. For an easily readable presentation of this classification of generic bifurcations, 
and of the corresponding normal forms used in semiclassical applications, we refer to the textbook 
of Ozorio de Almeida [16] . Bifurcations occurring in Hamiltonian systems with discrete symmetries 
have been investigated in [171 dHJ dSl [20] ; the TCB was, however, not mentioned in these papers. 
In [20] it has been shown that all other non-generic bifurcations occurring in such systems can be 
described by the generic normal forms given in [16], except for different bookkeeping of the number 
of fixed points which is connected not only to an m-tupling of the period, but also to degeneracies 
of the involved orbits due to the discrete symmetries. For the TCB this is not the case: it requires 
a normal form that is not in the generic list of [TB]. We derive an appropriate normal form for the 
TCB, starting from the general criteria given in [TT], and find it to correspond to that given in the 
literature for non-Hamiltonian systems [211 [22] . We use this normal form to develop the uniform 
approximation needed to include transcritically bifurcating orbits in the semiclassical trace formula. 
In a specific example that includes a TCB, we show numerically that our result allows to reproduce 
the coarse-grained quantum-mechanical density of states with a high precision. 

In the nonlinear and semiclassical physics community, there exists an occasional belief that 
non-generic bifurcations occur only in systems which exhibit discrete symmetries (time-reversal 
symmetry being the most frequently met in physical systems). The examples of TCBs which 
we present in this paper are obtained in a class of autonomous Hamiltonian systems with mixed 
dynamics; starting from the famous Henon-Heiles (HH) Hamiltonian [23] we change the coefficient 
of one of its cubic terms and add other terms destroying some or all of its discrete symmetries. All 
the TCBs that we have found involve one straight-line libration belonging to the shortest "period 
one" orbits. Our formal investigations therefore focus on the class of two-dimensional Hamiltonians 
containing a straight-line librational orbit. In this class of systems the TCB is, in fact, found to 
be the typical isochronous bifurcation of the librating orbit. The isochronous pitchfork bifurcation 
(PFB), however, which in Hamiltonian systems with time-reversal symmetry (such as the standard 
HH system) is the most frequently met non-generic bifurcation, is the exception here. We show how 
under a specific perturbation the PFB can unfold into a saddle-node bifurcation (SNB) followed 
by a TCB. In a specific example, we demonstrate that the TCB can exist in a system without 
any discrete (spatial or time-reversal) symmetry, thus proving that the above-mentioned belief is 
incorrect. 

Our paper is organized as follows. In Sec. [2] we compile results of Janich [111 [T2] relevant for 
our investigations. Starting from two-dimensional symplectic maps, we define a class of "crossing 
bifurcations" to which the TCB and the isochronous PFB belong. We discuss various criteria and 
properties of these bifurcations and give some useful formulae for the specific case of a bifurcating 
straight-line libration. The mathematically less interested reader may skip Sec. [2] and jump directly 
to Sec. [31 where we present numerical examples of the TCB and their characteristic features in 
the generalized Henon-Heiles Hamiltonians. We also study there various types of unfoldings of the 
TCB under perturbations of the Hamiltonian. In Sec. [H we discuss the semiclassical trace formula 
for the density of states of a quantum Hamiltonian, and present the uniform approximation by 
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which bifurcating periodic orbits can be included. In Sec. 14.41 we present a semiclassical calculation 
of the density of states in a situation where the TCB occurs between two of the shortest periodic 
orbits, and demonstrate the validity of the uniform approximation by comparison of the results with 
those of a fully quantum-mechanical calculation. In App.[6]we derive the appropriate normal forms 
for the TCB and the isochronous PFB which are needed in semiclassical applications. In App. [7] 
we briefly discuss the stability exchange of two orbits in a "false transcritical bifurcation" which 
actually consists of a pair of close- lying pitchfork bifurcations. 



2 Mathematical prerequisites 

In this section we present results of Janich jTH [T^ which are relevant for our investigations. We 
shall only quote theorems and other results; for readers interested in the mathematical proofs or 
other details, we refer to the explicit contents of [HI [12]. 



2.1 Poincare map and stability matrix 

We are investigating bifurcations of periodic orbits in two-dimensional Hamiltonian systems. They 
are most conveniently investigated and mathematically described by observing the fixed points on a 
suitably chosen projected Poincare surface of section (pss)Q Since the PSS here is two-dimensional, 
we describe it by a pair of canonical variables {q,p)- The time evolution of an orbit then corresponds 
to the two-dimensional Poincare map 

{q,p)^{Q,P), (2) 

where {q,p) is the initial and {Q, P) the final point on the PSS. Fixed points of this map, defined by 
Q = q, P = p, correspond to periodic orbits. We introduce e as a "bifurcation parameter" which in 
principle may be the conserved energy of the system or any potential parameter, normalized such 
that a bifurcation occurs at e = 0. Here we specialize to the energy variable by defining 

e = E-Eo, (3) 

where Eq is the energy at which the considered bifurcation takes place. We assume that the 
bifurcating orbit returns to the same point on the PSS after one map ([2]), so that Q = q, P = p; in 
this paper this will be called a "period one" orbit. We shall only study its isochronous bifurcations 
and hence only consider the non- iterated Poincare map. 

The map ([2]) is symplectic and thus area conserving in the {q,p) plane, and may be understood 
as a canonical transformation: 

Q = Q{q,p,e), P = P{q,p,e). (4) 

Janich [llj has given a classification of bifurcations of fixed points in two-dimensional symplectic 
maps, which we shall summarize in the following. We use his notation Q^, Pu for partial derivatives 
of the functions Q and P, respectively, with respect to u: 

O P - ^ (5) 



^ With "projected" we mean the fact that we ignore the value of the canonically conjugate variable (e.g., 
Py), to the variable (e.g., y) that has been fixed (e.g., by y = yo) to define the true mathematical PSS which 
lies in the energy shell. In the physics literature, it is standard to call its projection (with py = 0) the PSS. 
Due to energy conservation, the value of Py on the unprojected PSS can be calculated uniquely, up to its sign 
which usually is chosen to be positive, from the knowledge of q,p,yo and the energy E through the implicit 
equation E = H{q,yo,p,py), where H{x,y^Px,Py) is the Hamiltonian in Cartesian coordinates. 
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where u is any of the three variables q,p or e. Analogously Qqq, Pqpe, etc., denote second and higher 
partial derivatives. Due to the symplectic nature of the determinant of the first derivatives of 
Q and P with respect to q and p is unity: 

detf ^5,(9,^,^) QjiQ.P,^) )^^^ (g) 

We consider an isolated "period one" orbit with fixed point {q,p) = (0,0), for all values of e 
where it exists, and denote it as the A orbit. Its stability matrix is then given by 

M,(,)-( Qp(0,0,e) \ 

^^(^)-U,(0,0,e) P,(0,0,e) j- 

At e = 0, where the orbit undergoes an isochronous bifurcation, Ma(0) has two degenerate eigen- 
values +1, so that trMA(O) = 2. 

Henceforth we shall omit the arguments (0,0,0) in the partial derivatives of Q and P which - 
unless explicitly mentioned otherwise - will always be evaluated at the bifurcation point. When we 
need some of these partial derivatives at p = q = but at arbitrary values of e, we shall denote 
them by Qp{e) etc. When no argument is given, e = is assumed. We thus write 

The slope of the function tr MA(e) at e = (coming from a side where the orbit A exists) becomes, 
in this notation, 

trM^(O) =Q,e + Ppe. (9) 
By a rotation of the canonical coordinates q,p it is always possible to bring Ma(0) into the form: 

Ma(0)=(J ^l')' 

We shall henceforth assume that the coordinates have been chosen such that (llOp is trueH Then, 
with ^ one finds easily the determinant derivative formula [11] 

Qqu + Ppu = QpPqu (U = q,P, e) , (11) 

and dni) takes the simpler form 

trM:^(0) = QpP,,. (12) 

The total fixed point set 

F ■= {{q,P, e) I Q{q,P, e) = Q, P{q,P, e) = P} (13) 

is the inverse image of the origin (0, 0) in M2 under the map {Q — q, P — p) whose Jacobian matrix 
at (0, 0, 0) is 

Qq ~ 1 Qp Qe \ / Qp ^ (14) 



Pg Pp-l P, J V 



^In some cases one may find that Ma(0) has the transposed simple form in which Qp = and Pq ^ 0. In 
this case one may simply exchange the coordinates by a canonical rotation Q P, P ^ —Q (and q ^ p, 
p — > —q) in all formulae below and in App. 16.11 The case Qp — Qp = is exceptional and occurs only for 
harmonic potentials (cf. p4]). 
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In the generic case, 7^ and J has rank 2. This leads to the only generic isochronous bifurcation 
according to Meyer [TO], the saddle- node bifurcation (SNB) (also called "tangent bifurcation"). 
For this bifurcation, the fixed-point set F (1130 is a smooth one-dimensional manifold, consisting of 
two half-branches tangent to the q axis at the bifurcation point with slopes trM^(O) = itcxD. The 
orbit A then exists either only for e < or only for e > 0; no other orbit takes part in such a 
bifurcation. 

Following Janich [ITj, we speak of a rank 1 bifurcation, when the Jacobian J in (|14p has rank 
1, which is the case for 

= 0. (15) 

Then, after a suitable (e-dependent) translation of the p variable, J can always be brought into the 
form 

/ - 1 Qp \ / Qp 

V Pq Pp-1 J V 

We shall formulate all following developments in the suitably adapted coordinates {q,p), for 
which the form ()16p holds, and discuss only rank 1 bifurcations. 




2.2 Crossing bifurcations of isolated periodic orbits 

A rank 1 bifurcation for which the Hessian 

at (0,0,0) is regular and indefinite, i.e., for which detK = PqqP^^ — P^^ < 0, shall be called a crossing 
bifurcation. Janich showed [Tl] that a necessary and sufficient criterion for an orbit A to undergo 
a crossing bifurcation at e = is for the slope trM^(e = 0) to be finite and nonzero. With (jlOp and 
(fT2]) we see that 

Pge + (18) 

for crossing bifurcations. It follows that if the orbit A undergoes a crossing bifurcation at e = 0, 
it exists on both sides of a finite two-sided neighborhood of e = 0. Janich also showed that for 
such a bifurcation, the total fixed-point set F (|13p is the union ^4 U i? of two smooth 1-dimensional 
submanifolds intersecting at the bifurcation point. The set A is the set of fixed points (0, 0, e) of 
the A orbit; we shall call it the fixed-point branch A. The set B is the fixed-point set of a second 
orbit B which takes part in the crossing bifurcation. 

We shall discuss here only two types of crossing bifurcations: transcritical and fork- like bifur- 
cations. Their properties are specified in the following two subsections. A rank 1 bifurcation with 
a regular and definite Hessian K, i.e., with detK > 0, is sometimes called an "isola center" (cf. the 
normal form for the isola center in one-dimensional Hamiltonians at the end of Sec. 16. ip . Here the 
total fixed-point set F consists of the single isolated point (g,p, e) = (0,0,0). 
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2.2.1 Transcritical bifurcation (TCB) 

A transcritical bifurcation (TCB) occurs when, in the adapted coordinates {q,p) for which (jlOj) 
holds, one has 

P,,/0. (19) 

Then, there exists another isolated periodic orbit B on both sides of e = 0, forming a fixed-point 
branch B intersecting that of the orbit A at e = with a finite angle. The functions trMA(e) and 
trMB(e) have opposite slopes at the bifurcation: 

tr M^(0) = - tr M^3(0) . ( "TCB slope theorem") (20) 

In the scenario of a TCB, the orbits A and B simply exchange their stabilities and no new orbit 
appears (or no old orbit disappears) at the bifurcation. 

Note: Assume that the orbit A is a straight-line libration, chosen to lie on the y axis, so that 
the Poincare variables are q = x, p = px (see Sec. 12.31 below). Then, if the system is invariant 
under reflexion at the y axis, such a reflexion leads to P{q,p,e) = —P{—q,—p,e). Therefore, 
Pqq{q = 0,p = 0, e = 0) = Pgq = 0, and the bifurcation cannot be transcritical. The simplest 
possible crossing bifurcation then is fork-like (see next item). In short: Straight-line librations 
along symmetry axes cannot undergo transcritical bifurcations. 

2.2.2 Fork-like bifurcation (FLB) 

A fork-like bifurcation (FLB) occurs when one has 

Pqq = 0, P,<My^O. (21) 

Then, there exists another isolated periodic orbit B, either only for e > or only for e < 0. The 
fixed-point set of B consists of two half-branches intersecting the set A at e = at a right angle. In 
the adapted coordinates corresponding to p^ . one may parameterize the set B by {q,PBiq), ^siq)) 
and finds 

p's(0) = e's(0) = 0, e'ij(0)/0. (22) 

Although trMB(e) is not a proper function of e, a limiting slope trMg(O) ^ can be defined for 
both half-branches of the set B in the limit e ^ 0, coming from that side where they exist, and be 
shown [llj to fulfill the relation 

trM^(O) = -2trMA(0) . ( "FLB slope theorem") (23) 

In the same limit, the curvature of the set B at the bifurcation point is given by 

In the pertinent physics literature, this bifurcation is often called the (non-generic) isochronous 
pitchfork bifurcation (PFB). Note that here the two half-branches of the set B correspond to 
two different periodic orbits. They can be either locally degenerate (to first order in e), or globally 
degenerate due to a discrete symmetry (reflexion at a symmetry axis or time reversal). 

In the generic PFB corresponding to Meyer's classification [lUj . the fixed point scenario near 
e = is identical with that of the FLB. However, here the two fixed points of the set B correspond 
to one single orbit B which has twice the period of the primitive orbit A. In fact, the fixed-point 
branch A crossing the line trMA = 2 is that of the iterated Poincare map: the generic PFB is 
period doubling. The existence criterion (12 ip and the relations (1220 - ()24p for the B orbit hold 
here, too 1 251. 
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2.3 Some explicit formulae for straight-line librations 
2.3.1 Definition of the librational A orbit 

We now specialize to straight-line librational orbits in two-dimensional autonomous Hamiltonian 
systems, defined by Hamiltonian functions 

Ho{x,y,pa:,Py) = ^ (pI+pI) + V{x,y) (25) 

with a smooth potential V{x, y). Straight-line librations form the simplest type (and so far the only 
one known to us) of periodic orbits in Hamiltonian systems that undergo transcritical bifurcations. 
Let us choose the direction of the libration to be the y axis and call it the A orbit. The potential 
then must have the property 

1^ (»•='' = » <^«> 

for all y reached by the libration. The A orbit, which we assume to be bound at all energies, then 
has x{t) = x{t) =0 for all times t, and its y motion is given by the Newton equation 

dV 

m + ^{0,y{t))=0 y{t)=yA{t,e), (27) 

where yAit, e) is henceforth assumed to be a known periodic function of t with period ^^(e). For the 
A orbits in the (generalized) HH potentials discussed in the following section, the function y^(t,e) 
can be expressed in terms of a Jacobi-elliptic function [26]. We choose the time scale such that 
2/^(0, e) is maximum with 

yA(0,e)=0 Ve. (28) 

A suitable choice of Poincare variables is to use the surface of section defined by y = 0, and the 
projected PSS becomes the {x,px) plane, so that we define q = x, p = px- We again assume that 
the orbit A is isolated and exists in a finite interval of e around zero. The fixed-point branch A is 
thus again given by the straight line {qA,PA,G) = (0,0, e) in the {q,p,e) space. 

In [12j Janich has given an iterative scheme to calculate the partial derivatives Qq, Qp, etc. 
for this situation for any given (analytical) potential V{x, y) with the above properties. To this 
purpose, one has first to determine the fundamental systems of solutions ('^i,'^2) and {r]i,rj2) of the 
linearized equations of motion in the x and y directions, respectively: 

m + VxAO,yA{t,e))m = 0, (29) 
v{t)+Vyy{0,yA{t,e))ri{t) = 0, (30) 

with the initial conditions 

6(0) 6(0) \^( ??i(o) V2{0) A 1 
6(0) 6(0) J V ^i(o) m{0) J Vol 

For simplicity, we do not give the argument e of the 6(0 and r?i(t), but we should keep in mind that 
they are all functions of e. In (p9]) . ([30]) the subscripts on the function V{x,y) denote its second 
partial derivatives with respect to the corresponding coordinates. In the formulae given below, we 
denote by Vi{t), Vij{t), etc., with i,j G {x,y) the partial derivatives taken along the A orbit, i.e., 
at X = 0, y = yAit,e) as in ([29]) . ([30]) . evaluated at the bifurcation point e = 0. If the partial 
derivatives have no argument, they are taken at the period TA{eo), i.e., Vy = V^(T4(eo)) etc. 

Knowing the five functions yAit,e = 0) and 6(0) ^i(0 (^ = 1,2) at e = 0, all desired partial 
derivatives of Q{q,p,e) and P{q,p,e) at {q,p,e) = (0,0,0) can be obtained by (progressively re- 
peated) quadratures, i.e., by finite integrals over known expressions including these five functions, 
partial derivatives of V{x,y), and the functions obtained at earlier steps of the scheme (whereby 
the progression comes from increasing degrees of the desired partial derivatives). 



Ve. (31) 
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2.3.2 Stability matrix of the A orbit 

We note that the equation (f29l) is nothing but the stabihty equation of the A orbit, since the by 
definition are small variations transverse to the orbit. In the standard literature, ()29p is also called 
the "Hill equation" (cf., e.g., [14^ I27j). The stability matrix Ma at the bifurcation of the A orbit is 
therefore simply given by 

with Ta = Ta{€ = 0). Its eigenvalues must be Ai = A2 = +1, as seen directly from (fTO]l . The 
solutions £,i{t,e) of (I29p are in general not periodic. But at the bifurcations of the A orbit, where 
trM^ = +2, one of the ^j(t,e = 0) is always periodic with period Ta (or an integer multiple m 
thereof) [27j and describes, up to a normalization constant depending on e, the transverse motion 
of the bifurcated orbit at an infinitesimal distance e from the bifurcation (cf. [Ml [Ml [2S])- 

2.3.3 Slope of the function trMA(e) at e = 

Here we give the explicit formulae, obtained from [12j, for the slope tr M^(0) — Qqe ~1~ Ppe^ see (|9]) . 
of the function trMA(e) at the bifurcation. The quantities Qq^ and Pp^ are given, in terms of the 
potential V{x^ y) in ([25]) and the other ingredients defined above, by 

Qi- = t^pMTa)-^ I' \xxy{t)[QpUt)-QqUt)]Ut)m{t)dt, (33) 

\^y) Jo 

= ^-^^QpmiTA)-^ [ \..yit)[Ppm-P<iUt)]Ut)viit)dt. (34) 

[^y) ^y Jo 

In the adapted coordinates where trMA(O) has the form (jlOp with Qq = Pp = 1 and Pq = 0, the 
slope becomes 

tr M:,(0) = Qp Pq, = ^^-^ Qp m{TA) -^Qp V,,y{t) e?(t) rji{t) dt . (35) 
K^y) ^y Jo 

For the case that trMA(O) has the transposed tridiagonal form with Qp = and Pq ^ 0, the formula 
becomes 

\t)rjiit)dt. (36) 



tvM'j,{0)=PqQp,= -L^Pqm{TA) + ^Pq [ \..y{t) 
[Vy) Vy Jo 



An independent derivation of ()33p - ()36p is given in [30j, where it is shown that the first terms 
are due to the variation of the A orbit's period Ta with e, whereas the integral terms are due to the 
e dependence of the functions (t) . 



2.3.4 Criterion for the TCB 

For a bifurcation to be transcritical, we need Pqq ^ 0. From |12j we find the following explicit 
formula for Pqq 

Pqq = - [ ^V,,,{t)f^{t)dt, (37) 

Jo 

which also yields explicitly the parameter b in its normal form given in ()96|) below. 

If the potential is symmetric about the y axis, then Vxxx{t) is identically zero and the TCB 
cannot occur, as already stated in Sec. 12.2.11 above. However, even if Vxxxit) is not zero, special 
symmetries of the function ^i(t), in combination with that of Vxxxi't)^ can make the integral in (j37p 
vanish. An example of this is discussed in Sec. 13.3.51 
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3 TCBs in the generalized Henon-Heiles potential 



3.1 The generalized Henon-Heiles potential 

For our numerical studies, we have investigated the fohowing family of generalized Henon-Heiles 
(GHH) Hamiltonians: 



(38) 



Here a is the control parameter that regulates the nonlinearity of the system, and 7, /3 are parameters 
that define various members of the family. The standard HH potential [23] corresponds to 7 = 1, 
13 = 0. It has three types of discrete symmetries: (i) rotations about 27r/3 and 47r/3, and (ii) 
reflections at three corresponding symmetry lines, which together define the Cst, symmetry, and 
(iii) time-reversal symmetry. There exist three saddles at the critical energy E* = l/6a^ so that 
the system is unbound and a particle can escape if its energy is ii^ > E* . For 7 7^ 1, /? 7^ 0, the 
spatial symmetries are in general broken (except for particular values of 7 and /?) and only the 
time-reversal symmetry is left. There still exist three saddles, but in general they lie at different 
energies. There is always a stable minimum at x = y = 0. 

It is convenient to scale away the nonlinearity parameter a in (j38|) by introducing scaled variables 
x' = ax, y' = ay and a scaled energy e = E/E* = 6a^E. Then (|38|) becomes 

h = e = Qa^E = 3 {pi, + p^, + x'"^ + y'^) _ 2 y'^ + 6 (7 x'^y' + (3 y'^x') , (39) 

so that one has to vary one parameter less to discuss the classical dynamics. (For the standard HH 
potential with 7 = 1, /3 = 0, the scaled energy e then is the only parameter.) For simplicity, we 
omit henceforth the primes of the scaled variables x', y'. 

Before we discuss the periodic orbits in the system (j39p . let us briefly recall the situation in the 
standard HH system in which all three saddles lie at the scaled energy e = 1. 



3.1.1 Periodic orbits in the standard HH potential 

The periodic orbits of the standard HH system have been studied in [28} [31} [32], and their use 
in connection with semiclassical trace formulae in [33l [Ml [35l [36l [SZl [38]. We also refer to f23] 
(section 5.6.4) for a short introduction into this system, which represents a paradigm of a mixed 
Hamiltonian system covering the transition from integrability (e = 0) to near-chaos (e > 1). 

In Fig. [T] we show the trace tr M of the stability matrix M, henceforth called "stability trace" , 
of the shortest orbits as a function of e. Up to energy e ~ 0.97, there exist [3Tj only three types 
of "period one" orbits [in the sense defined after ([3|)]: 1) straight-line librations A along the three 
symmetry axes, oscillating towards the saddles; 2) curved librations B which intersect the symmetry 
lines at right angles and are hyperbolically unstable at all energies; and 3) rotational orbits C in the 
two time-reversed versions which are stable up to e ~ 0.89 and then become inverse-hyperbolically 
unstable. While the B and C orbits exist at all energies, the orbits A cease to exist at the critical 
saddle energy e = 1 where their period becomes infinite. 

When |trM| > 2 or < 2, an orbit is unstable or stable, respectively. When trM = 2 it either 
undergoes a bifurcation if the orbit is isolated, or it belongs to a family of degenerate orbits in the 
presence of a continuous symmetry. The latter is seen to occur in the limit e — > 0, where the orbits 
A, B and C all converge to the family of orbits of the isotropic two-dimensional harmonic oscillator 
with U(2) symmetry. The A orbits undergo an infinite sequence of (non-generic) isochronous PFBs, 
starting at e ~ 0.97 and cumulating at e = 1. At these bifurcations an alternating sequence of 
rotational orbits (labeled R) and librational orbits (labeled L) are born. This bifurcation cascade. 
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the R and L type orbits, and their self-similarity have been discussed extensively in [26[I28|. In Fig. [T] 
and in the text below, we indicate their Maslov indices (needed in semiclassical trace formulae, see 
Sec. HI by suffixes to their labels, which allows for unique bookkeeping of all orbits. (Only the first 
two representatives R5 and Lg of the orbits born along the bifurcation cascade are shown in Fig. [T] 
by the dashed lines.) At each bifurcation, the orbit A increases its Maslov index (which is 5 up to 
the first bifurcation) by one unit. Only the first three bifurcations can be seen in the figure; the 
others are all compressed into a tiny interval below e = 1. As has been observed numerically in 
[Ml [23 [32]) trM^ becomes a periodic function of the period Ta in the limit e — > 1. [It can actually 
be rigorously shown that, asymptotically, trM^(T^) — > —2.68042 sm{^/3 Ta) in this limit |3Uj .] 




I 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

e 

Figure 1: Trace of the stability matrix M of the "period one" orbits in the standard HH potential, 
plotted versus scaled energy e. The suffixes indicate their Maslov indices. Only the first two (R5 
and Lg) of the orbits born at an infinite sequence of isochronous PFBs of the A orbit, cumulating 
at the saddle energy e = 1, are shown (dashed lines). 

As is characteristic of isochronous PFBs (cf. Sec. 12.2.2]) . the new-born orbits come in degenerate 
pairs due to the discrete symmetries: the two librational L orbits are mapped onto each other under 
reflection at the axis containing the A orbit, and the two rotational R orbits are connected by time 
reversal. Note that although the A orbit ceases to exist for e > 1, all R and L type orbits bifurcated 
from it exist at all energies e > 1. For some new orbits appearing there, we refer to the literature 
[321 138j : in this paper we shall not be concerned with them. 

3.1.2 Periodic orbits in the generalized HH potential 

For 7/0, /? / 0, there exist in general three different saddles at scaled energies eo, ei and 62, and 
three different straight-line periodic orbits, labeled A, A' and A", oscillating towards the saddles. 
In general, there are three curved librational orbits B, B' and B" (not necessarily existing at all 
energies) intersecting the three A type orbits at right angles, and there is always a time-reversally 
degenerate pair of rotational orbits C going around the origin. It is rather easy to see that the three 
A type orbits always intersect each other at the minimum of the potential located at the origin 
{x,y) = (0,0). The equations of motion for the Hamiltonian ()38p in the Newton form are (in the 
scaled variables corresponding to a = 1) 
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x + x{l + 2jy)+ Py"^ = 0, 



For a straight-line orbit librating through the origin we have y = ax which, inserted into (140 p . yields 
a cubic equation for the slope a: 

/3a3 + (27 + l)o2-2/3a-7 = 0. (41) 

For the rest of this paper, we limit the parameters to the range 7 > and /? > 0. Then, ()4ip has 
always real roots that are in general different. In the right panel of Fig. [2] below, we have shown 
the six shortest ( "period one" ) librations obtained numerically for 7 = 0.6, (3 = 0.07, including the 
three straight-line orbits A, A', A" intersecting at the origin. 

Further analytical analysis is cumbersome except for the following special cases: 
/J = 0, 7 = 1 (standard HH) : 

Two of the slopes are ai^2 = ±1/V3; the third is oo = 00 corresponding to the orbit along the y axis 
with x{t) = 0. The three saddles lie at {x,y) = (0,1), (-\/3/2, -1/2), and (\/3/2, -1/2), forming 
an equilateral triangle with side length \/3 ; its sides (and their extensions) form the equipotential 
lines for e = 1. The periodic orbits are those discussed in Sec. I3.1.T1 
/3 = 0,7/1: 

The rotational C^^ symmetry is broken, but the reflection symmetry at the y axis is kept. Corre- 
spondingly, we find two degenerate orbits A', A" with opposite slopes 01^2 = ±\/7/(27 + 1). There 
is a horizontal equipotential line at yi = 1/2 = — 1/27 with scaled energy ei = 62 = (3 + l/7)/472 
that contains two saddle points symmetrically positioned at xi^2 = ±\/(2 + I/7) /27. At low ener- 
gies, there is only one B type orbit intersecting the y axis at a right angle; two further orbits B' 
and B" appear through bifurcations at higher energies (see examples in Sec. 13. 2p . For 7 > there 
is a third A orbit librating along the y axis (ag = 00) towards a third saddle at (0, 1) with energy 
eo = 1. The equipotential line for e = ei^2 consists of the horizontal line at yi^2 = "1/^7 and two 
branches of a hyperbola. For 7 > 1 the hyperbola branches lie symmetrically about the y axis, each 
intersecting the horizontal line at one of the two symmetric saddle points. For < 7 < 1, they lie 
symmetrically about a horizontal line at y* = (1 + 37)/47, the lower of them intersecting the line 
y = yi,2 at the two symmetric saddle points. 

The limiting case /3 = 0, 7 = yields a separable and hence integrable system with only one 
saddle at (0, 1) at energy eo = 1 and one A orbit (with a = 00). We do not discuss this system here, 
but refer to [37] in which it is investigated both classically and semiclassically in full detail. 

3.2 Examples of transcritical bifurcations and their properties 

As mentioned above, we have restricted the parameters 7 and /? in the GHH potential ()38p to be 
positive (or /? = 0). We find that, depending on the values of /? and 7, at least one or two of the 
straight-line orbits A, A', or A" can undergo a TCB with a partner of the curved librational orbits 
B, B', or B". In the following, we shall first show two examples and then discuss characteristic 
properties of the TCB. In Sec. 13.31 we shall study its stability and its unfoldings. Some of the 
numerical results can easily be understood analytically in terms of the normal forms of the various 
bifurcations and their unfoldings. These are discussed in detail in App. [6] and shall be referred to 
in the following text. 
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3.2.1 Two examples 



As a numerical example, we choose 7 = 0.6, (3 = 0.07. The three saddle energies are eo = 0.993 for 
the A orbit, ei = 2.81 for the A' orbit, and 63 = 3.74 for the A" orbit. In the left panel of Fig. [2] we 
show the stability traces trM(e) of the shortest orbits. In the right panel we display the shapes of 
these orbits in the (x, y) plane. The orbits B' and B" are created in a SNB at Cth ~ 1.533 and do not 
exist below this energy; at high energies they are hyperbolically unstable with increasing Lyapunov 
exponents. Contrary to the standard HH system, only the A orbit is stable at low energies, while 
the orbits A' and A" leave the e = limit unstable and cross the critical line trM = +2 at some 
finite energies etch and ej^^ to become stable. At higher energies, all three A type orbits undergo an 
infinite PFB cascade as in Fig. [H each of them converging at its saddle energy. (We do not show 
here the R and L type orbits born at these bifurcations.) 

It is between the pairs of orbits A', B and A", B' that we here observe two TCBs. They occur 
at the energy etcb = 0.854447 between the orbits A' and B, and at e'^^^ = 1.644 between the orbits 
A" and B'. The situation near e'^^^^ actually displays an example of a slightly broken PFB which will 
be discussed in Sec. 13.3.51 




Figure 2: Left panel: Stability traces of A and B type orbits in the GHH potential with 7 = 0.6, 
(3 = 0.07, plotted versus scaled energy e. The three saddles are at eo = 0.993, ei = 2.81, and 
62 = 3.74. Right panel: Shortest orbits, projected on the (x, y) plane. Orbit A is evaluated at 
e = 0.99 just below its saddle (eo = 0.993); all other orbits are taken at e = 2. The line types 
correspond to those in the left panel. 
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Figure 3: TCB of a degenerate pair of orbits A', A" and B', B" at etch = 0.4889 in the GHH 
potential with 7 = 0.75, /3 = 0. The degeneracy is due to the reflection symmetry about the y axis. 
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Another example of a TCB is shown in Fig. O obtained for the GHH potential with 7 = 0.75 
and P = 0. This potential is symmetric about the y axis and therefore the pairs of orbits A', A" 
and B', B" are degenerate, lying opposite to each other with respect to the y axis. The crossing 
happens at etcb = 0.4889 and exhibits the same features as those discussed in the first example. 

3.2.2 Characteristic properties of the TCB 

We now discuss some of the properties of a TCB and compare our numerical results to their 
analytical predictions from the normal form of the TCB. For this purpose, we take the example at 
etcb = 0.854447 seen in Fig. where the orbits A' and B bifurcate transcritically. Their crossing 
is shown in Fig. U] on an enlarged scale in the upper left panel, where the numerical results for 
trM(e) are displayed by crosses (orbit A') and circles (orbit B). We see that the graphs of trM(e) 
cross the critical line trM = 2 with opposite slopes. Their Maslov indices, differing by one unit, are 
exchanged at the bifurcation (see Sees. HI and 14. 2p . The upper right panel displays the numerical 
action difference AS = 83 — Sa' (circles), where the action of each periodic orbit (po) is, as usual, 
given by 

V = ^P-dq. (42) 

In the lower panels, we show the shapes of the orbits in the (x,y) plane below (left) and above 
(right) the TCB. The B orbit is seen to have passed through the A' orbit at the bifurcation. The 
lengths of both orbits increase with energy e. 




-0.5 0.0 0.5 -0.5 0.0 0.5 

X X 



Figure 4: TCB in the GHH potential with 7 = 0.6, /? = 0.07. Orbits A' and B exchange their 
stabilities (and Maslov indices a = 4,5) at etct = 0.854447. Upper left: trM versus energy e; crosses 
(A') and circles (B) are numerical results, solid lines the local prediction P5|) . Upper right: action 
difference AS* versus e; circles are numerical results and the solid line the local prediction (j99p . 
Lower panels: shapes of the crossing orbits in the (x, y) plane before and after the bifurcation. 
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The normal form of the TCB is derived and discussed in Sec. 16.2.31 From it, one can derive the 
local behavior of the actions, periods, and stability traces of the two orbits in the neighborhood of 
a TCB. For small deviations e = c (e — etcb) (with c > 0) from the bifmxation energy, the stability 
traces go like trM(e) = 2 it 2ae, and the action difference of the two orbits like AS'(e) = — 
(see App. 16.2.31 for the meaning of the other parameters). These local predictions, given in Fig. U] 
by the solid lines, can be seen to be well followed by the numerical results. 

The crossing of the graphs tr M(e) of the two orbits at the bifurcation energy etcb with opposite 
slopes is a characteristic feature of the TCB (see Sec. 12.2.1]) . Since the fixed points of the two orbits 
coincide at the bifurcation point, their shapes must be identical there. In the present example, the 
orbit B is a curved libration; the sign of its curvature is changed at the bifurcation, as illustrated 
in the two lower panels of Fig. [H 

We note that a completely different mechanism of stability exchange of two orbits, which happens 
through two close-lying PFBs, has been described in [39]. The stability diagram may then appear 
like that in the upper left of Fig. [H if the crossing point is not analyzed with sufficient numerical 
resolution. Such a "false transcritical bifurcation" will be briefly discussed and illustrated in the 
appendix of our paper. 



3.3 Stability and unfoldings of the TCB 

Since the TCB is not a generic bifurcation according to Meyer's list [10], we now address the 
question under which circumstances it can exist and what its structural stability is. The GHH 
systems discussed here have time-reversal symmetry, and it is therefore of interest to study the 
stability of the TCB under perturbations of the Hamiltonians that destroy this symmetry. In 
this context, it is important to note that a detailed mathematical study [IT], in which all generic 
bifurcations in systems with time-reversal symmetry are classified, does not mention the TCB; the 
same holds also for |20J. So far we have only found TCBs which involve a straight-line libration. On 
the basis of the results presented in Sec. [21 we believe that in the class of all Hamiltonian systems 
containing a straight-line librating orbit, the TCB is actually the generic isochronous bifurcation 
of the librating orbit. Therefore, if we find a perturbation of the GHH system that destroys the 
time-reversal symmetry but preserves a straight-line libration, the TCB should also exist there. 
This will be demonstrated in Sec. 13.3.41 for a specific example. 

A general Hamiltonian H{x, y,Px,Py) supports the existence of a straight-line libration - which, 
without loss of generality, may be chosen to lie on the y axis - if the following conditions are fulfilled: 

dH dH 

— (0,y,0,p,)=0, _(0,y,0,p,)=0. (43) 

In the following we will first show how some TCBs are destroyed under perturbations that violate 
the conditions (j43p . and how they unfold. We find two types of unfoldings which are also discussed 
in [SI [21] for TCBs in non-Hamiltonian systems. In the first scenario, the TCB breaks up into 
SNBs. In the second scenario, no bifurcation is left in the presence of the perturbation and the 
functions trM(e) approach the critical line trM = 2 without reaching it, so that one may speak of 
an avoided bifurcation. These scenarios can be described by the extended normal forms given 
in Sec. 16.2.41 We then also investigate perturbations that fulfill the criteria ()43p . allowing for the 
existence of TCBs in systems with or without any discrete symmetries. 
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3.3.1 Addition of a homogeneous transverse magnetic field 



We first discuss the addition of a homogeneous magnetic field B = b^Bq to the Hamiltonian ()38p 
which is transverse to the (x, y) plane of motion. This is a situation that is frequently set up in 
experimental physics and gives us one important way of breaking of the time-reversal symmetry. 
The momenta pi (i = x,y) in (j38p are replaced by the standard "minimal coupling", 

Pi ^ Pi--Ai, A = ^ (r X B) , (44) 
c 2 

where A is the vector potential and e the charge of the particle. This adds the following perturbation 
to the Hamiltonian: 

SH{x,y,px,Py) = ^{xpy-yp:^) + ^(^^^ (x^ + y^) , (45) 

which breaks the time-reversal symmetry of (j38p due to the linear terms in px and py, but also 
breaks the straight-line libration condition (|43p . 

As an example, we choose the GHH potential with 7 = 0.5, /3 = 0.1. Here the saddle energy 
for the A' orbit is ei = 3.83; the other saddles are at eo = 0.9852 and 62 = 6.35. In Fig. [5] we show 
the stability traces trM(e) of the orbits A' and B' with and without magnetic field. For Bq = 
(triangles and dashed-dotted lines), these orbits A' and B cross at euf = 1.42665 in a TCB like in 
the examples discussed above. For Bq (circles and solid lines), they rearrange themselves into 
pairs A'4/B5 and B4/A5 colliding in SNBs according to the prediction (jl04p of the normal form 
(jl03p . in which n is taken proportional to the value of Bq. 




Figure 5: Unfolding of a TCB by a transverse magnetic field in the GHH potential with 7 = 0.5, 
P = 0.1. Shown is trM versus scaled energy e. Dashed lines and triangles: prediction (|98p and 
numerical results for field strength i?o = of the unperturbed TCB; solid lines and circles: local 
prediction (I104p (with an adjusted value of k) and numerical results for Bq = 0.0001. 
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3.3.2 Destruction of the TCB by a perturbation of the potential 

Another example of the same unfolding of a destroyed TCB is shown in Fig. [6l Here the unperturbed 
GHH potential is the same as that used in Fig. [3] above, which is symmetric about the y axis. This 
time we apply a perturbation of the potential alone 

5V{x,y) = Kx'y'^ , (46) 

whereby x\ y' are rotated Cartesian coordinates such that the bifurcating A' orbit lies on the 
y' axis. Clearly, this perturbation does not fulfill the conditions ()43|) (expressed in the rotated 
coordinates) and in fact destroys the original TCB of the orbits A and B' shown in Fig. [3l the 
same fate happens also to the pair A" and B" of orbits. We see in Fig. [6] that, again, the original 
pairs of orbits on either side of the unperturbed TCB rearrange themselves such as to destroy each 
other in two pairs of SNBs, each according to the prediction (jl04p of the corresponding normal 
form. Since the effective perturbation strengths are different in the two original directions of the 
A' and A" orbits, the splitting between the two pairs of SNBs is slightly different. A problem 
arises with the nomenclature of the perturbed orbits, which is somewhat ad hoc, since all perturbed 
orbits have become rotations. In the square brackets in the figure we indicate the names of the 
unperturbed orbits, of which A', A" are straight-line and B', B" curved librations (their stability 
traces are shown in Fig. [3] above) . The stability traces of the perturbed orbits change drastically 
at the original bifurcations, but approach those of the unperturbed orbits sufficiently far from the 
bifurcations. The insert in the upper left of Fig. [6] illustrates one possible unfolding of a destroyed 
isochronous PFB (that seen at e = 0.34 between the orbits B and B'/B" in Fig. and will be 
commented in Sec. 13.3.51 below. 




0.0 0.1 0.2 0.3 0.4 0.5 

e 



Figure 6: Unfolding of the TCB shown in Fig. [3] under the perturbation ([^6]) with k = 0.0001 (see 
text for details). The labels i n brac kets [] correspond to the orbits of the unperturbed system in 
Fig. El For the insert, see Sec. 13.3.51 



16 



3.3.3 An avoided TCB 



In Fig. [7] we give an example of an avoided bifurcation. We start again from the same example as 
in Fig. [3l but now we apply the following perturbation: 

5H{x,y,px,Py) = k' x'^pyi , (47) 

again in the same rotated coordinates as for the perturbation (|46|) above. By construction, this 
perturbation does fulfill the libration-conservation conditions (|43|) . expressed in the rotated coordi- 
nates {x',y') for the orbit A', so that the TCB of the orbits A' and B' survives. It will be discussed 
in more detail in the next subsection. The perturbation (I47p destroys, however, the TCB of the 
original orbit pair A" and B" at etcb = 0.489 and is seen to lead to an avoided bifurcation of the 
perturbed orbits, which are now called A" and B and shown by the heavy dashed lines. Their 
stability traces follow the local behavior (jl06p predicted by the normal form (jlOSp . Again, our 
nomenclature for the new orbits is not strict; the perturbed B" orbit has, for e > 0.489, become a 
portion of the new orbit B. As in Fig. [6l the graphs trM(e) of the perturbed new orbits approach 
the unperturbed ones far from the bifurcations. 
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Fig ure 7: Unfolding of the TCB of the orbit pair A" and B" shown in Fig. [3]under the perturbation 
(j47p with k! = 0.01 into an avoided bifurcation of the new orbit pair A" and B near e = 0.489 
(shown by the heavy dashed lines). The surviving TCB of the orbit pair A' and B' (solid and thin 
dashed line, respectively) is commented in Sec. 13.3.41 



3.3.4 TCB in a system without any discrete symmetry 

We now come to our last, and perhaps most interesting, example: a TCB in a system without any 
discrete symmetry. It is shown in Fig. [7] by the solid line for the orbit A' and the thin dashed line 
for the orbit B'. It is the same as the TCB shown in Fig. [3] after applying the perturbation ()47|) 
that has been explicitly constructed so as to preserve the straight-line libration condition ()43p in 
the rotated coordinates x', y' . Here y' is the direction of the A' orbit. Thus, the libration A' in the 
perturbed system is identical to that in the unperturbed GHH potential (7 = 0.75, /? = 0). The 
orbit B', however, which in the unperturbed GHH system is a curved libration similar to that shown 
in Fig. SI has now become a rotation except at the TCB point. While it was originally created, 
together with its symmetry-degenerate partner B", in an isochronous PFB at e = 0.34 from the 
original B orbit (see Fig. [3]), this PFB is destroyed under the perturbation (I47p . and the perturbed 
B' orbit is now created at a SNB at e = 0.343. Its stable lower branch is that which crosses the 
unchanged A' orbit transcritically at the slightly shifted new bifurcation energy etch = 0.4886. 

The shapes of this perturbed B' orbit in the rotated {x',y') plane are shown in Fig. [8l on the 
left side in the energy region between its creation at e = 0.343017 and its TCB at e = 0.4886 where 
it is stable, and on the right side for the energies e > 0.4886 where it is unstable. Its librational 
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shape at e = 0.4886, where it is identical to the A' orbit, is shown in both panels of the figure (note 
their difi'erent scales!). 




Figure 8: Shapes of the B' orbit, shown by the dashed line in Fig. [71 in the rotated {x\ y') plane at 
different energies. Left panel: stable region, at the energies e = 0.343 (creation in SNB) - 0.4886 
(TCB point). Right panel: unstable region, at the energies e = 0.4886 (TCB point) - 4.0. Note 
the different scales in the two panels. 




Figure 9: Same as Fig. [H]but in the rotated momentum space {Px'iPy')- 

In Fig. [9] we present the shapes of the B' orbit in the rotated momentum space {Px',Py')- Here 
the orbit appears as a figure-8 type rotation, except at the TCB where it must be a straight-line 
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libration, as the A' orbit, also in momentum space. It should be noted that qualitatively, the shapes 
in momentum space are the same for all transcritically bifurcating B type orbits discussed in this 
paper, even if they remain curved librations in coordinate space. 

One may interpret the perturbation (I47p as the first-order expansion of a weak inhomogeneous 
magnetic field with strength proportional to k'. In a homogeneous field, for which the full pertur- 
bation is given in (I45p . the Lorentz force tends to turn all straight-line orbits into curved librations. 
In the system perturbed by (j47p . the Lorentz force of this inhomogeneous magnetic field is canceled, 
at least to lowest order, by the geometry of the total potential which tends to curve the librations 
the other way round. We leave it to the interested reader to speculate whether this scenario finds 
applications in accelerator physics, where one may want to produce straight-line trajectories in an 
inhomogeneous magnetic field. 

3.3.5 Creation of a TCB in the unfolding of a PFB 

In this section we will show how a TCB can be created by perturbing a PFB and how its existence 
may depend on particular symmetries. Two characteristic unfoldings of isochronous PFBs in one- 
dimensional dynamical systems have been described in |8lj9j, which correspond to the "universal 
unfolding" of the isochronous PFB discussed extensively in |21| . We find the same unfoldings for the 
isochronous PFBs of the straight-line librations in the (G)HH potentials; one of them is of particular 
interest here as it leads to a TCB. The corresponding normal forms are given in App. 16.2.61 

In the first scenario, the original parent orbit does not change its stability, thus avoiding the 
bifurcation, and a pair of new orbits is created at a SNB. One of these new orbits takes the role 
of the original parent orbit after the bifurcation, and the perturbed parent orbit takes the role of 
one of the new orbits created at the original PFB. Examples of this scenario can be seen in Fig. [6] 
(inserted close-up) and in Fig. [3 as results of two different perturbations of the same original PFB 
seen in Fig. [3]at e ~ 0.34. 

The second scenario is the unfolding into a SNB of a pair of new orbits, followed by a TCB 
of one of these orbits with the original parent orbit. An example of this has already been pointed 
out in Fig. [2] (left panel) to occur near e ~ 1.62, where the orbits B' and B" bifurcate from the A" 
orbit. In the following we shall further illustrate this unfolding by explicitly perturbing a PFB in 
the standard HH system. 

We start from the HH system, i.e., (i38]l with 7 = 1, /? = 0, and add the following perturbation 
to the potential 

5V{x,y) = ^5x\ (48) 

which destroys both the symmetry and the reflection symmetry at the y axis. It therefore 
affects the cascade of isochronous PFBs of the linear A orbit along the y axis (cf . Sec I3.1.ip . The 
perturbation (j48p is chosen such as to preserve the straight-line libration condition (j43p . so that 
the A orbit still exists in its presence. To ensure the presence of a TCB in the perturbed system, 
we must fulfill the condition Pqq 7^ given in (jl9p . An explicit expression for the quantity Pqq in 
terms of the (total, perturbed) potential V is given in (fHTj) . [In the integrand of ([57]) . the function 
yxxx{x,y) is taken along the A orbit with x{t) = 0, y = yA{t); see Sec. 12.31 for details and notation.] 
Since Vxxx becomes nonzero with the perturbation (I48p . the occurrence of a TCB is possible. But 
Vxxx 7^ is not sufficient to ensure Pqq 7^ 0: this will also depend on the symmetry of the function 
^i(t) appearing in the integrand of the quantity Pqq in (f37|) . 

Now, as discussed in [26\ [28], the functions (,i{t) describe the x motion (transverse to the A 
orbit) of the new orbits created at the successive PFBs of the A orbit. These functions are periodic 
Lame functions with well-known symmetry properties. As it turns out, ^i(t) of the L type orbits 
born at every second PFB of the cascade are even functions of t, where t = is the time at which 
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yAit) is maximum; whereas those of the R type orbits born at every other bifmxation are odd. The 
result is that Pqq becomes zero at the R type bifurcations, in spite of Vxxx / 0, while Pgg 7^ for the 
L type bifurcations. Consequently, it is only at the L type bifurcation energies that a TCB can exist 
in the perturbed system. Our numerical investigations have confirmed that under the perturbation 
()48p all R type bifurcations remain, indeed, unbroken PFBs with unchanged stability traces trM(e) 
to first order in 5, while the L type bifurcations are broken up as discussed above. 

In Fig. [To] we show the creation of the orbits Lg and L'g from the A orbit in the HH system 
under the perturbation (08]). In the unperturbed HH system, these orbits are born as a degenerate 
pair from a PFB at the energy eg = 0.986709235 (cf. as also seen in Fig. [TJ Here the PFB has 

been broken according to the second scenario described above, unfolding into a TCB of and 
Ly/g at precisely the same critical energy etcb = eg, and a SNB at eg — Aeg ~ 0.98703 where L7 and 
L'g are born. The thin dash-dotted line gives the slope of tr M(e) of the L orbit at the TCB which is 
minus that of the A orbit, as is characteristic of a TCB. The thin dotted line gives the slope of the 
original degenerate pair Lg, L'g created in the unperturbed HH system; this slope is minus twice 
that of the parent A orbit, as is typical of a PFB [see (j23p in Sec. 12.2.2] . The same scenario is found 
at all successive L type bifurcations. In the perturbed system, the orbit pairs L, L' are no longer 
degenerate, since the reflection symmetry at the y axis is broken by (j48p . The graphs trM(e) seen 
in this figure are predicted by the normal form (|113p . as discussed explicitly in App. 16.2.61 
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Figure 10: Creation of orbits Lg, L'g in a broken PFB in the HH potential under the perturbation 
(I48p with 6 = 0.5. The TCB of the orbits A and L occurs at the energy eg = 0.986709235 of the 
PFB of the degenerate pair L, L' in the unperturbed HH potential [2B]. (See text for the thin lines.) 
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4 Inclusion of a TCB in the semiclassical trace formula 



4.1 The semiclassical trace formula 

Our investigations have been largely motivated by the use of periodic orbits in the semiclassical 
description of the density of states g{E) in a quantum system with discrete spectrum {Ei} 

g{E)=Y,5{E-Ei). (49) 

i 

Initiated by Gutzwiller (see [13j and earlier references therein), the periodic orbit theory (POT) 
(see also |15t [2^ W\\ for general introductions) states that the oscillating part of the quantum 
density of states is given, to leading order in h, by a semiclassical trace formula of the form 

5g{E) = Ao{E) cos - | a,,] . (50) 

The sum goes over all periodic orbits {po) of the classical system (including their repetitions), Spo{E) 
are their actions (fi2|) and dpo their Maslov indices (see [I3l Il2l HSj ISj and [2Sj, appendix D, for 
details). The Apo{E) are semiclassical amplitudes which depend on the nature of the orbits. For 
orbits that are isolated in phase space, the amplitudes were given by Gutzwiller p!3] as 



Apo{E) = \ , ^"'"^^^ , (51) 

^\/|det[M^„(^)-I]| 



in terms of their primitive periods Tppo{E) = dSppo{E)/dE and stability matrices Mpo{E); here I 
is the unit matrix with the same dimension as Mpo. For systems with continuous symmetries, in 
which most periodic orbits come in degenerate families (in particular in integrable systems), explicit 
expressions for the Apo{E) have been derived by various authors [i5] . 

One problem with the Gutzwiller trace formula in mixed systems, where stable and unstable 
periodic orbits coexist, is the divergence of the amplitudes ()5ip occurring at bifurcations. Remedy is 
given by uniform approximations, introduced by Ozorio de Almeida and Hannay [l6] (see also [16J) 
and further developed by several authors both for codimension-one [25 1 147 ^ ^8] and codimension-two 
bifurcation scenarios |37l I49| . In the following sections, we will discuss the uniform approximations 
and derive its appropriate form for a pair of transcritically bifurcating orbits. 

The trace formula ()50p does not converge in mixed systems and most chaotic systems, in which 
the number of periodic orbits proliferates exponentially with increasing length, so that the sum- 
mation over all orbits typically cannot be performed (see [14^ 140] ) . In our study, we coarse-grain 
the density of states by convolution with a normalized Gaussian with width AE, so that only the 
shortest orbits with periods Tpo ^ h/AE contribute to the sum [231 [501 EH- Although the finer 
details of the spectral information hereby are averaged out, the coarse-grained density of states 

1 

^/nAE ^ 



9ae{E) = ^^ye-^^-^^^y^^^' (52) 



still exhibits its gross-shell structure, provided that AE is not chosen too large. The correspondingly 
coarse-grained trace formula reads [24J 

5g^E{E) = Y: ApoiE) e-[^- ^^/^^l^ cos f _ 1 ^ J , (53) 



h 2 

po ^ 



where it can be seen that the additional exponential factor suppresses the contribution of longer 
orbits. This version of the POT has found many applications to gross-shell effects in finite fermion 
systems (see [231 [52] for examples) . 
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4.2 Uniform approximation for bifurcating orbits 

In this section we sketch the derivation for the combined contribution of a pair of bifurcating orbits 
A and B to the semiclassical trace formula (|50p for the density of states. Since the individual 
amplitudes in the form (j5ip given by Gutzwiller diverge at the bifurcation, one has to go one step 
back in their evaluation and transform the trace integral to the phase space [25\ [53l . After doing 
the integration along the primitive A orbilH with action Sa{E), the remaining part of the trace 
integral is over the Poincare surface of section in the variables Q and p transverse to the A orbit: 

Sg{E) = 3f?e e^^^^^^-'^''^ J dQ J dp C{Q,p, e) e^^^^'P'^l (54) 

The action function in the phase of the integrand is given by 

S{Q,p,e) = S{Q,p,e)-SA{e)-Qp, (55) 

where S{Q,p,e) is the generating function of the canonical transformation 1^ that describes 
the Poincare map, and 5'A(e) is the action integral (j42|) of the A orbit as a function of the control 
parameter e. By virtue of the canonical relations of the generating function 

the stationary condition of the function S in the {Q,p) plane for any fixed e, 

dS dS 

^(Qo,]5o,e) = — (Qo,]5o,e) = 0, (57) 

yields Pq = po and Qq = qq, so that the stationary points {Qo,Po,e) = {qo,Po,€^) of the phase 
function (I55p are the fixed-point branches of the map and hence correspond to the periodic orbits. 
[Note that, by construction, S{qo,po, e) = along the fixed-point branch of the A orbit.] 

The amplitude function C{Q,p,€) in (j54p is given ^25j in terms of the generating function 
S{Q,p,e) by 

1 (99 

W,P,e) = ^^(Q,^',e)- (58) 

Note that 

g^{Q,P, e) =: T{Q,p, e) = Ta{E) + ^(Q,P, e) , (59) 

where Ta{E) = -^Sa{E) is the period of the A orbit. 

In principle, the integration over Q and p in (I54p is limited to that domain of the {Q,p) plane 
which is accessible under energy conservation. However, in the spirit of the stationary-phase ap- 
proximation (including its extensions below) we expect that, due to the rapidly oscillating phase 
of ()54p in the semiclassical limit S ^ h, the main contributions to the integral come from small 
regions around the stationary points of the function S{Q,p,€). Assuming that the fixed points of 
A and the other orbit(s) taking part in the bifurcation are situated in the interior of this domain, 
and that no other bifurcations happen nearby, we may extend the integrals over both Q and p from 
— oo to -|-oo. 

Sufficiently far away from the bifurcation point e = 0, so that the orbits A and B are isolated, 
the stationary-phase integration of ()54p will yield precisely the contributions of the isolated orbits 
A and B to the standard Gutzwiller trace formula ()50p with their individual amplitudes (I5ip . Near 



^recall that we only consider primitive "period one" orbits and their isochronous bifurcations here. 
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the bifurcation, the stationary-phase approximation fails and one has to include higher than second- 
order terms in Q and p in the function S{Q,p,e). The simplest solution [16j is to use a truncated 
Taylor expansion of S{Q,p,e) in all three variables, keeping only the minimum number of terms 
necessary to be able to reproduce locally the fixed-point branches of all the orbits taking part in a 
given bifurcation. These truncated forms of S{Q,p,e) are the normal forms which are discussed 
in App. m 



4.3 Uniform approximation for the TCB 

Equipped with the normal form of the TCB given in App. 16.2.31 we now calculate the contribution 
of a pair of periodic orbits A and B undergoing a TCB to the semiclassical trace formula. We 
follow closely the treatment of [2^, where uniform approximations for the generic bifurcations 
corresponding to [10\ [T6] were derived. 

Since (j54p is invariant under canonical transformation {Q,p) iQ',p'), we may think of the 
variables Q,p to be the adapted coordinates for which (|85|) and the equations given thereafter are 
valid. We are therefore allowed to insert for S{Q,p, e) the normal form derived in Sec. 16.2.51 for the 
TCB, in order to derive the uniform approximation to the trace formula which includes the orbits 
taking part in the TCB. 

We will do this in two steps. First, we evaluate ()54p only at e = 0. This yields the so-called local 
uniform approximation in the spirit of Ozorio de Almeida and Hannay [46j. In the second step, 
we use the full normal form ([96]) and the corresponding functions C{Q,p, e) defined by ([58]) to find, 
after some suitable transformations, the global uniform approximation in the spirit of |25l 1471 H8] . 
The latter yields asymptotically the Gutzwiller trace formula for the orbits A and B sufficiently far 
from the bifurcation. 

We now use for S{Q, e) the normal form (llOip of the TCB (omitting the tilde on the variables 
Q,p,€ and on b). Using the relation ([3]), the function T in (I59p becomes 

f{Q,e) = TA{E)-Q\ (60) 

After the elementary p integration, yielding a complete Fresnel integral, we obtain for ([54 



dg{E) = ^5Reei^^(^)-*f("^+i)[T4(^)F(6,e)-G(6,e)], (61) 

7r?iv 27r/i 

where we have defined the two following one-dimensional integrals 

/oo 
dQe-fA^Q'+'Q'), (62) 
-oo 

/oo o 
dQQ^e-i^''^'+''^'^ = ih^F{b,e). (63) 
-oo 

Using the substitution Q = x — Q with Q = e/3b, we obtain 

F{b,e) = 27T(J^Y\TA^^'^Ai{-z'), z'= ^ , (64) 



^3|6|y ' " (3|6|)4/3;i2/3 ' 

where Ai is the Airy function (see [57], 10.4) and 



3 



AS(.) = -|^. (65) 
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Using the r.h.s. of (j63|) to calculate G{b, e), we finally obtain for the level density 



V2 



^Wl^ (3|6|)V3 



^.(^) + iJ)Ai(-.') + .^S'"Ai'(-.') 



(66) 



4.3.1 Local uniform approximation 



We first give the result ()66p for e = 0, using the known value [57] of Ai(0), to find the local uniform 
level density at the bifurcation energy £^0^ 



Ta{Eo)T{1) 
7rV6^nV6|6|i/3 



cos 



lsA{Eo)-laA-l 



(67) 



which contains the combined contribution of both orbits A and B taking part in the transcritical 
bifurcation. An explicit expression for the calculation of the normal form parameter b is given in 
([37|) of Sec. 12.31 The result ([67j) looks identical to that obtained in [25] for the generic SNB. The 
reason is that the normal form for this bifurcation is [25J S{Q,p,e) = —eQ — bQ^ — (7p'^/2 which 
for e = gives, of course, the same result as the normal form (|1U1|) . Note that the power 7/6 of ^ 
in the denominator is by 1/6 higher than in the semiclassical amplitude ()5ip of an isolated orbit. 



4.3.2 Global uniform approximation 

The result ()67p gives the correct semiclassical amplitude of the bifurcating pair of orbits A and 
B only locally at the bifurcation, i.e., for e = 0. We want, however, to know it also away from 
the bifurcation, and in particular, also in the limit where it can be written as a sum of the two 
individual contributions of the isolated orbits A and B to the standard Gutzwiller trace formula 
(j5Up . To achieve this, we note that if we use the asymptotic forms of the Airy function and its 
derivative in (j66p for \z'\ ^ 1, we obtain two terms that formally look like contributions to (j5Up 
with amplitudes of the form (I5ip . but with the actions Spo{E), periods Tpo{E), and stability traces 
tr Mpo{E) replaced by their expansion to lowest order in e, as found from the normal form and given 
in dM]) and (l99]) . 

The next intuitively obvious step is therefore to rewrite the asymptotic form of (I66p in terms of 
the locally expanded quantities Spoie), Tpo{e), and trMpo(e) of the two orbits (po = A,B) and then 
to replace them by the correct functions Spo{E), Tpo{E), and ticMpo{E) found numerically for the 
isolated orbits away from the bifurcation. This step has been rigorously justified in [25] by some 
appropriate transformations and need not be repeated here. The calculation goes exactly like that 
presented in [23] for the case of the SNB on that side where both orbits are real. The reason is that 
although the normal forms of the two bifurcations are different, they lead to identical integrals after 
a translation in the integration variable Q. 

The result is the following uniform contribution of the two bifurcating orbits to the Gutzwiller 
trace formula: 

= (m _ _ M£) 3i„ (M _ Ai'(-4 . (68) 
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The quantities occurring in (|68p are defined as 



z = (3^/2)2/3, e = 2-15^-5^1, 

5 = ^ (•S'a + Sb) , CT = ^ (CTA + (Tb) , 

A = ^{Aa + Ab), AA = ^{AA-AB)sign{SA-SB), (69) 

all to be taken at the energy E, where Spo{E) and cTpo are the actions and Maslov indices, respec- 
tively, of the isolated periodic orbits on either side of the bifurcation, and Apo{E) > are their 
Gutzwiller amplitudes (f5T]) . 

At the bifurcation (e = z = ^ = 0), the result ()68p reduces to the local uniform approximation 
()67p . Far enough away from the bifurcation, it goes over to the contribution of the isolated orbits 
A and B to the standard Gutzwiller trace formula. Indeed, expressing the Airy function in terms 
of Bessel functions as [57] 

Ai(-^) = lV^[Jl/3{0+J-l/3{0], 

Ai'(-^) = ^^[J2/3(0- ^-2/3(0], (70) 



and using their asymptotic form 




^ cos fe - J - t) for ^ » 1 ' (71) 



vr^ V" 2 4 



we obtain from ([68]) for ^ ^ 1, i.e., for \Sa — Sb\ ^ 2/i, the sum of the isolated Gutzwiller 
contributions to the trace formula 

6g{E)= J2 AoiE)cos(^^-^apX (72) 

po=A,B ^ ^ 

with the amplitudes Apo{E) given in (f^Tj) . 

For the reason given above, the result (|68p looks identical to that given in for the SNB on 
that side where the two orbits are real. The present result holds on both sides of the TCB and can 
easily be seen to yield asymptotically the result (f72]) on both sides, with the roles of the orbits A 
and B and their Maslov indices properly exchanged. 



4.4 Numerical test 

We present here a numerical calculation of the density of states, both quantum-mechanical and 
semiclassical, of the GHH Hamiltonian (j38p with 7 = 0.6, /3 = 0.07, whose shortest orbits and 
stability traces are shown in Fig. [2j The parameter a in the Hamiltonian (j38p has been chosen as 
a = 0.04. The three saddles then lie at the energies Eq ~ 103, Ei ~ 293 and E2 — 390 in units such 
that the spacing of the harmonic-oscillator spectrum reached in the limit E ^ equals huj = 1. 
(See Fig. [2] for the scaled energies eo, ei and 62 of the saddles.) 

The spectrum of the quantum-mechanical Hamiltonian corresponding to (I38p has been obtained 
by diagonalization in a two-dimensional harmonic oscillator basis. Strictly speaking, the spectrum 
is not discrete since the system has no lower bound. However, for energies below the three saddles, 
the tunneling probabilities are exponentially small. In principle, the semiclassical trace formula 
can also be applied in the continuum region above the saddles, if the (complex) energies of the 
resonances are used to calculate the density of states. For a detailed discussion of this situation, 
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we refer to [HH] where a semiclassical calculation has been successfully performed for the standard 
HH system up to twice the saddle energy. In the present system, the discrete energies obtained 
for E'o ~ -E' ~ 150 in the numerical diagonalization turn out to be good approximations to the real 
parts of the resonance energies, while the imaginary parts of the resonances are still negligible. 

For our present test, we have chosen a coarse-graining width in (I52p of = 0.6. This allows 
us to restrict the summation over the periodic orbits (po) to the primitive ("period one") orbits; 
including second or higher repetitions does not affect the numerical results within the resolution of 
the lines presented in the figure below. As we can see in Fig. [2] (left panel), there exist only five 
"period one" orbits in the system below the scaled energy e ~ 1.5 corresponding to ~ 156. 
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Figure 11: Oscillating part 5g{E) of the density of states in the GHH potential (f38l) with 7 = 0.6, 
P = 0.07 and a = 0.04, plotted versus the unsealed energy E. The quantum-mechanical (qm) 
results are shown by the solid lines (identical in all three panels), different semiclassical (scl) 
approximations are shown by the crosses. Both qm and scl results have been coarse-grained by a 
Gaussian with width AE = 0.6, see (|52p . Top: the five shortest (primitive) orbits A, B, C, A' and 
A" are included by the standard Gutzwiller trace formula ()5ip . ()53p for isolated orbits. Note the 
divergence near E = 89, corresponding to etcb = 0.854447 where the TCB of orbits A' and B occurs 
(cf. Fig. U]). Middle: same as in the top panel, but the two crossing orbits A' and B are omitted. 
Bottom: same as in the middle panel, but now the crossing orbits A' and B are added in the global 
uniform approximation (|68p . 

In Fig. \TT\ we show the oscillating part of the level density obtained for a = 0.04 as a function 
of the unsealed energy E, up to = 135 which corresponds to e ~ 1.3. The solid lines display the 
coarse-grained quantum-mechanical result which is the same in all three panels. In order to extract 
the oscillating part of the quantum density of states (f52]) , we subtracted its Strutinsky averaged part 
which here corresponds to the Thomas-Fermi approximation [23]. The crosses, connected by dashed 
lines, represent the semiclassical result in various approximations. The regular fast oscillations with 
period ~ 1 on the energy scale E come from the common average action Spo{E) of the leading 
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periodic orbits, which becomes the action Sho{E) = 271 E of the harmonic oscillator in the limit 
S — > 0. The beat-like slow variation in the amplitude of 6g{E) is due to the interferences of the 
periodic orbits and can be captured by the semiclassical trace formula ()53p . 

In the top panel of Fig. \TT\ the five orbits A, B, C, A' and A" are included in the trace formula 
(I53p with their Gutzwiller amplitudes (I5ip . Although they qualitatively reproduce the main trends 
of the beating density of states, they overestimate it. One can clearly see the divergence at £" ~ 89 
corresponding to the scaled energy etcb = 0.854447, where the TCB of the orbits A' and B occurs. 
(The divergences due to the PFB sequence of the orbit A near e = 0.993 <^ E = 103.5 cannot be 
seen with this resolution.) The center panel shows the same semiclassical result, but omitting the 
contributions of the bifurcating orbits A' and B. Clearly, the agreement with quantum mechanics 
is not good even far from the bifurcation, showing that these orbits always play a role. In the 
bottom panel, the orbits A, C and A" are again included as isolated orbits with the amplitudes 
()5ip . while the combined contribution of the bifurcating orbits B and A' is included in the global 
uniform approximation given in (I68p . The agreement between semiclassics and quantum mechanics 
is now excellent, demonstrating the adequacy of the uniform approximation. The fact that the 
isolated-orbit approximation in the top panel does not work even far away from the bifurcation 
shows that the orbits A' and B do not become isolated enough in the energy region shown; i.e., the 
asymptotic form (|72p of the uniform approximation is not reached. This could already be expected 
from the fact that the stability traces trM(e) of these orbits stay very close to +2 for all e < 1.2, 
as can be seen in Fig. [2j 

In the energy limit — > (not shown in Fig. [TTI) , the present semiclassical approximations are 
not appropriate due to the integrable limit of the harmonic oscillator. A corresponding uniform 
approximation for the standard HH potential has been derived in [36]. It can be generalized in a 
straightforward manner to the GHH systems, following the lines of [36], but this would lead beyond 
the scope of the present paper. 
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5 Summary 



We have discussed transcritical bifurcations (TCBs) of periodic orbits in non-integrable two-dimen- 
sional autonomous Hamiltonian systems. We have first discussed the mathematical aspects of the 
TCB, making use of recent studies by Janich [11^ I12j . We then have, with the help of numeri- 
cal examples in generalized Henon-Heiles (HH) systems, discussed their phenomenology and their 
unfoldings under perturbations. We have shown, in particular, that a TCB may also exist in a 
system without any discrete symmetry, although it does not belong to Meyer's list [TO] of generic 
bifurcations. The reason is our restriction to systems containing straight-line librations. In such 
systems, the TCB appears to be the generic isochronous bifurcation of the straight-line libration 
while its isochronous pitchfork bifurcation (PFB) represents an exception expressed by the condi- 
tion Pqq = in (I19p . Using this condition and the explicit expression (I37p for the calculation of 
Pqq for Hamiltonians of the form (I25p , we have exploited a special type of unfolding of the PFB to 
construct a perturbation of the standard HH system under which the TCB occurs. 

So far, we have only encountered TCBs of straight-line librational orbits. In most examples, 
the second orbit that takes part in the TCB with the straight-line libration is also a librational 
orbit, though not along a straight line. That this need not be so has been shown in the example in 
Fig. [HI where the second orbit is a rotation (except, of course, at the TCB where it coincides with 
the straight-line libration) . This system, in the presence of the momentum-dependent perturbation 
()47|) , is the only example of a Hamiltonian which does not have the form psp and in which we found 
a transcritically bifurcating straight-line orbit. From the general criteria given in Sec. 12.21 however, 
we see no a priori reason why rotational orbits should not undergo TCBs as well. Furthermore, it 
is obvious that a given straight-line orbit can always transformed into a more complicated one by 
suitable canonical transformations. The inverse question - if an arbitrary nonlinear periodic orbit 
can be canonically transformed into a straight-line libration which can bifurcate transcritically - 
might also have a positive answer, but we see no way of proving or testing this, nor can we give a 
nontrivial example, since periodic orbits (except straight-line librations) in non-integrable systems 
usually cannot be given analytically. 

Finally, we have constructed a global uniform approximation for the inclusion of transcritically 
bifurcating periodic orbits in the semiclassical trace formula for the quantum density of states. A 
numerical comparison with the fully quantum-mechanical calculation of the coarse-grained density 
of states yields excellent agreement. The normal forms of the TCB and the isochronous PFB have 
been derived in the appendix, where we also point to a "false TCB" which is the result of a stability 
exchange between different orbits via an intermediate periodic orbit through a pair of PFBs. 
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6 Appendix 1: Normal forms 



In singularity theory (see, e.g., [21]) and catastrophe theory (see, e.g., [55]) it is standard to classi- 
fy bifurcations according to their normal forms. In App. 16.11 we briefly discuss normal forms for 
isochronous bifurcations in non-Hamiltonian one-dimensional systems and give the explicit forms 
for the TCB, the PFB and the saddle-node bifurcation (SNB). 

While non-Hamiltonian fields can always be transformed to normal forms by suitable coordinate 
transformations, the situation is more complicated for Hamiltonian fields [56]. Here the normal 
forms depend on pairs of canonical variables, like the Poincare variables {q,p) used in Sec. [2l and 
must be derived, for a given Hamiltonian and a given type of bifurcating orbit, from the generating 
function S{Q,p,e) of the Poincare map or, equivalently, from the function S{Q,p,e) given in (f^^ . 
by suitable canonical transformations. The strength of the normal forms - if they can be found - 
is that they are unique for each generic type of bifurcation and do not depend on the particular 
form of the Hamiltonian or the bifurcating orbit. However, the reduction of S{Q,p,e) to one of 
these generic normal forms is, according to Arnold [56], "generally not possible, and formal series 
for canonical transformations reducing a system to normal form generally diverge" . 

Nevertheless, for the generic bifurcations occurring in two-dimensional symplectic maps, as 
analyzed and classified by Meyer [lOj, normal forms suitable for semiclassical applications have 
been given in [THl US] 



6.1 Normal forms for one-dimensional non-Hamiltonian systems 

We follow here the book of Golubitsky and Schaeffer [2T] and use their notation. In a simple one- 
dimensional problem with a "state variable" x and a "bifurcation parameter" A, one may study the 
set of values (x. A) satisfying the equation 

5(x,A) = 0, (73) 

where g{x, A) is a differentiable scalar function of both arguments. Bifurcations of this set occur at 
critical points (xo,Ao) where 



= . (74) 

xo,Xo 



The TCB can be specified by the following criteria: at the critical point (xo,Ao), the function g 
must fulfill 

9 = 9x = 9x = 0, 9xx 7^0, det (fg = g^x g\\ - glx<^, (75) 

where the subscripts denote partial derivatives with respect to the corresponding variables. For the 
fixed points Xn = x^+i = x of the quadratic map we get the function 9{x,r) = rx{l — x) — x 
which fulfills the criteria (I75p with A = r at the critical point {xq = 0, tq = 1), so that a transcritical 
bifurcation must occur there. 

Normal forms are the simplest functions, usually taken to be polynomial forms in x and A, which 
obey the criteria for a given bifurcation. They can often be found by Taylor expansion of a given 
equivalence class of functions around the critical points, keeping the lowest necessary number of 
terms required to fulfill the given criteria. A valid normal form for the TCB is [22] 

gtcbix, A) = -Ax, (76) 



^ We point out a misprint that occurred in both [16] and [46]: the normal form for the generic SNB was 
erroneously given analogous to that in ([M)) below, which is the normal form of the TCB. For the SNB the 
first term should correctly be —eQ rather than —eQ^. 
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which fulfills the criteria (j75|) at (xq, Aq) = (0, 0). Normal forms are not unique. The following form 
is easily seen to be strongly equivalent (in the sense of [21]) to (j76]) : 



gsb{xA) = ^x' -X', (77) 

since it also fulfills the criteria ([75]) at (xq, Aq) = (0,0). Although the bifurcation corresponding to 
(j77p in [2T] is referred to as "simple bifurcation", it is identical to what we here call the TCBj 
The criteria for the isochronous PFB are 



9 = gx = gx = gxx = , s^a 7^ o , g^xx + , (78) 

and its standard normal form is 

gp/fe(3;,A) = - Ax. (79) 

For completeness we give here also the criteria for the SNB (or tangent) bifurcation (in [21j 
called "limit point"): 

5 = fx = , 5a 7^ , ffxx 7^ , (80) 

and its normal form: 

gsnh{x, A) = - A . (81) 
Golubitzky and Schaeffer [2^ also list the "isola center" bifurcation whose criteria are: 

g = gx = g\, gxxy^O, det cPg>0, (82) 

with the normal form 

gUx,X) = ±ix^ + X^). (83) 

In two-dimensional Hamiltonian systems, the isola center is, according to Janich's classification [llj . 
a rank 1 bifurcation for which the Hessian matrix K in ()17p is regular and definite. 



6.2 Normal forms for crossing bifurcations 

In order to find normal forms for the two types of crossing bifurcations discussed in this paper, which 
we need in the semiclassical trace formula (|54p . we follow the heuristic approach of determining a 
truncated Taylor expansion of S{Q,p, e) with the minimum number of terms necessary to describe 
the required properties of these bifurcations. To this purpose, we first establish relations between 
the partial derivatives of the functions Q{q,p, e) and P{q,p, e) in ([!]), and the partial derivatives of 
the function S{Q,p, e) for which we use the same notation as in 12.11 Translating the criteria given 
in Sec. l2.2l for the crossing bifurcations in terms of the partial derivatives of S, we can determine the 
normal forms of the TCB and the FLB. Although the formal transformations needed to arrive at 
these normal forms are not necessarily canonical, their use in the semiclassical trace formula can be 
justified by the fact that possible missing terms of higher order do not affect the results to leading 
order in h (cf. [ISl [23 [13 [IH] ) • 



^Note that in ^T], the name "transcritical" is used for a whole class of bifurcations, differently from our 
restricted use of the term. 
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6.2.1 Relations between Q,P and S and their partial derivatives 

From (I55p and ()56p we obtain the following basic relations 

= q -Sp{Q,p,e) , 
P{q,p,e) = p + SQ{Q,p,e). (84) 

We now take partial derivatives of these relations with respect to the variables q, p and e, in order 
to formulate the conditions for the crossing bifurcations discussed in Sec. 12.21 in terms of partial 
derivatives of the function S{Q,p,e). This procedure is simplified by the following step. The 
splitting lemma of catastrophe theory (see ^55j, pp 95 and 103) states that after a suitable (but 
perhaps not canonical) coordinate transformation, the function S{Q,p,e) can be split up in the 
following way 

S{Q,p,e) = S{Q,e)-^p^ a^O, (85) 

where S{Q,e) does not depend on p any more. In the suitably adapted coordinates q,p and Q,P 
for which (1850 is truell we obtain the following relations 

Qg{e) = l, Qp[e) = -Spp{e) = a , 

= SQQie) , Pp{e) =l + a Sggie) , (86) 

and tr M(e) becomes 

trM(e) = 2 + cT5QQ(e), (87) 

which is valid along the fixed-point branches of both orbits A and B. We also give some of the 
higher partial derivatives of Q and P at e = (valid in the adapted coordinates): 

Qqq = ; Qqp = ; ^qq = ^QQQ > Pqp = ^ ^QQQ 

and 

Pqe = SqQe , Pqqq - 



6.2.2 Criteria for the two crossing bifurcations 

We can now express the criteria for the two types of crossing bifurcations introduced in Sec. 12.21 
directly in terms of the parameter a and the partial derivatives of the function S{Q,€) defined in 
()85p . To have a bifurcation of the A orbit at e = 0, we must have 

Sqq{Q = 0, e = 0) = . {bifurcation of A orbit) (90) 

For the occurrence of a rank 1 bifurcation, we have the criterion (see the end of Sec. 12. ip 

P^ = SQe = . (rank 1 bifurcation) (91) 

Since {P, Q) = (p, q) = (0, 0) then is the fixed-point branch of the A orbit for all e, the function 
S{Q,e) must fulfill, due to ([8^ . the condition 

5^(0,6) =0 Ve. {fixed-point branch of A orbit) (92) 

The criterion for the occurrence of a crossing bifurcation is that the slope trM^(O) given by 
([12]) be nonzero, see also (fT8|) . We therefore need 

Spp = —cr / , SgQe 7^ . [crossing bifurcation) (93) 
^and for which Ma(0) has the form PH)) [see also the footnote after PU)) ] 
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The criterion for this bifurcation to be transcritical is 

Sqqq 7^ 0. [transcritical bifurcation) (94) 

For the occurrence of a fork-hke bifurcation, we must have 

Sqqq = , Sqqqq / . {fork-like bifurcation) (95) 

We are now ready to construct the simplest normal forms for the function S{Q,p,e), split like 
in dM}, that fulfin all the above criteria ([90]) - ([93]) and either or ([95]) . 

6.2.3 Normal form of the TCB 

For the transcritical bifurcation, the normal form obtained in this way is 

SiQ,p,e) = -eQ^-bQ^-^p\ (6/0) (96) 

with 6 = — I Pqq. An explicit formula for calculating Pqq and hence the parameter b is given in ()37p . 
The normal form (j96p corresponds to that of the TCB in non-Hamiltonian systems given in (j76p in 
App. [6m if we choose Sq{Q,p=0, e) = gtcbiQ,^), but it has, to our knowledge, not been discussed 
in connection with bifurcations of periodic orbits in Hamiltonian systems. 
The fixed-point branch of the B orbit is easily found to be 

2 % 
PB{e) = 0; QB{e) = --e ^ eB{Q) = --Q. (97) 

The stability traces of the two orbits are then found from (j87p to be 

trMA(e) = 2-2o-e, 

trMB(e) = 2 + 2o-e, (98) 

fulfilling the "TCB slope theorem" ([20|) . Along the branch B, the function S{Qb-,Pb-,() yields a 
contribution to the action of the B orbit. Noting that the contribution to the A orbit is zero, this 
yields the action difference of the two orbits: 

S{Qb,Pb, e)=AS = SB-SA = -^. (99) 

Note that a sign change of either a or e in (|96p simply corresponds to exchanging the orbits A 
and B, whereas a sign change of b does not affect the local predictions (f98l) and (f99]) . 

In the applications of the normal forms for semiclassical uniform approximations, one usually 
assumes a = ±1 (see, e.g., |47[ 125)). However, when starting from an arbitrary Hamiltonian, this is 
not automatically fulfilled. In fact, one sees directly from (I86p that a = Qp which a priori is not of 
modulus unity. But we can easily absorb the absolute value of cr by a canonical stretching (shear) 
transformation: {Q,p) (Q , p) specified by 

Q = Q/VW\, p = pVW\- (100) 

The normal form (j96p then becomes 

S{Q,p,e) = -~eQ^-bQ'-^p\ a = ±1 (101) 

with 

e = \a\e, b = \a\^/^ b . (102) 
In Sec. 14.21 we shall use the form (jlOip but omit the tilde on all variables and constants. 
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6.2.4 Normal forms for two unfoldings of the TCB 

We have found two scenarios for the destruction of a TCB by a perturbation K5H(x,y,px,Py) of 
the Hamiltonian, where k is a real parameter. In the first scenario, the bifurcation unfolds into a 
pair of SNBs lying opposite to each other on either side of the unperturbed bifurcation point e. This 
scenario can be described by adding to the normal form (j96p a term linear in Q with a negative 
sign: 

SiQ,p,€) = -fi^Q- eQ^-bQ^-^p^ (6 > 0) (103) 

It predicts the following local behavior of the stability traces: 

tr Ma,b = 2 ± 2cr\/e2 - 36^2 . (104) 

Between the two SNBs, which occur at e = ±V3bK, there are no real periodic orbits. For e < — VsEk 
and for e > V^k, the pairs of original orbits A and B join and "destroy" each other in the SNBs. 
Examples for this scenario are given in Sec. 13.3.11 where k is proportional to the strength of a 
homogeneous external magnetic field, and in Sec. 13.3.21 

In the second scenario, no bifurcation is left in the presence of the perturbation. The two pairs 
of orbits approach the critical line trM = 2 from both sides, come closest to it at the original 
bifurcation point e = 0, and then diverge from it again. We will call this an "avoided bifurcation". 
It can be described by a normal form identical to (|1U3|) . except for an opposite sign of the linear 
term: 

S{Q,p,€) = +K^Q- eQ^-bQ^-^p'^. (6 > 0) (105) 
This form predicts for the local behavior of the stability traces 

trMA,B = 2 ± 2cr\/e2 + 36^2 , (106) 
corresponding to an avoided bifurcation. An example of this unfolding of the TCB is given in 

Both types of unfoldings have been found also in non-Hamiltonian one-dimensional systems (see, 
e.g., [8]). According to [21], one may describe a "universal unfolding" of the TCB by the normal 
form 

SiQ,p,e) = 5Q-eQ^-bQ''-^p^ (6 > 0) (107) 

This corresponds to (|lU3p for 6 < and to p05|) for 6 > 0. The above two normal forms are, 
however, mathematically different in that they predict only one type of unfolding for both signs of 
the parameter k in each case, whereas ()107p would predict different unfoldings for the two signs 
of one and the same parameter 5. We emphasize that all three cases describe different (physical) 
phenomena; we have not found any realization of the type pOTp in our numerical studies of TCBs. 

6.2.5 Normal form of the FLB 

For the fork-like bifurcation (FLB), we arrive at the normal form 

S{Q,p,e) = -eQ'-aQ'-^p^ (a / 0) (108) 

with a = —-^Pqqq- This form is identical to that of the generic period-doubling PFB [TB] and 
corresponds, with SQ{Q,p=0,e) = gpfbiQ,^), to that of the isochronous PFB in non-Hamiltonian 
systems given in (f79]l . 

The two fixed-point branches of the B orbits here are given locally by 

PB{e)=0; Qeie) = ±V-e/2a ^ eeiQ) = -2aQ\ (109) 
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where the rightmost relation fulfills the conditions e^(0) = 0, e^(0) 7^ given in (|22p . The stability 
traces of the A and B orbits (on that side of e = where the latter exist) are found locally to be 

trMA(e) = 2-2c7e, 

trMB(e) = 2 + 4cTe, (110) 
fulfilling the "FLB slope theorem" (I23p . and their action difference becomes 

S{QB,PB,e)=^S = SB-SA = ^. (Ill) 

4a 

The same local behaviors (IllOh and (lllip hold also for the generic PFB |25j . 

Note that the B branches describing the bifurcated new orbits B only exist on that side of the 
bifurcation where e/a < 0. Changing the sign of e has the same effect as changing the sign of a. 
For the GHH potentials, all isochronous PFBs of the A type orbits have negative values of a. 

Changing the sign of a "mirrors" the bifurcation scenario at the line trM = +2, i.e., the 
stabilities of all orbits are exchanged from stable to unstable and vice versa. 



6.2.6 Normal form for unfoldings of the FLB 

We have found two scenarios for unfoldings of the FLB, which are known also to occur in non- 
Hamiltonian systems [SI [21]. Since the FLB locally is identical to the isochronous PFB, one may 
use its universal unfolding in one-dimensional non-Hamiltonian systems |21] as a guide and choose 
the normal form 

5(g,p,e) = -eQ2_aQ4^5Q + KQ3_|^2^ _^ ^-^^2) 

thus adding a linear and a cubic term in Q to (jlOSp . One type of unfolding is obtained for 6^0. 
Hereby a pair of orbits is created in a SNB on one side of a critical value ei of e (depending on 5 
and k), while a third orbit is present on both sides of e and does not undergo any bifurcation. This 
is the most usual unfolding of the FLB (and persistent in the sense of |21J)- Examples are given in 
Fig. [6] (inserted close-up) and in Fig. [71 

The second scenario, more interesting for our present investigations, is obtained from (|112p as 
the special case 6 = 0. Hereby the FLB (or isochronous PFB) is broken up into a TCB, occurring at 
the original bifurcation point e = 0, associated with a near-lying SNB. Examples of this unfolding 
are shown in Fig. [2] (left panel) and in Fig. [TOl The fixed-point set in this unfolding is found from 
Sp = p = and 

5Q(Q,e) = -2eQ - 4aQ^ - 3kQ2 = 0, (113) 

yielding the fixed-point branch {Qa,Qb) = (0,0) of the original A orbit, which is not affected by 
the perturbation with n ^ 0, and a set of new points {Qb,Pb = 0), where Qb are solutions of the 
equation 

4aQ^ + 3kQ + 2e = 0. (114) 

This corresponds to a parabola in the {Q,e) plane. Calculating the stability traces trMyi(e) and 
tiMsie), one obtains the graphs seen in Fig. [TUl The A branch is locally linear around e = 0. The 
B branches, realized here by the Lg, L7 and Lg orbits, form a tilted parabola with the following 
properties, (i) The symmetry line of the parabola has the slope trM^(O) of the original B orbits 
created at the unperturbed FLB, obeying the "FLB slope theorem" tr M^(0) = -2tr M^(0) given in 
(I23p . (u) The parabola intersects the A orbit in the new TCB with the slope — trM^(O) according 
to the "TCB slope theorem" ([20]). 



34 



7 Appendix 2: False transcritical bifurcations 



We briefly address here a mechanism of stabihty exchange which has been described in [39]. Hereby a 
pair of periodic orbits exchange their stabihties via an intermediate periodic orbit which is exchanged 
between the two other orbits through two oppositely oriented isochronous PFBs. On a moderate 
scale of the control parameter, the intermediate orbit may not be observed numerically, and the 
stability traces of the two main orbits appear to cross the critical line trM = 2 exactly like in a 
TCB. These are, however, "false transcritical bifurcations". 

We illustrate this in Fig. [121 calculated for the coupled quartic oscillator potential V{x,y) = 
(x^ + y^)/4 + ax^y^/2 which is non-integrable for all values of a except a = 0, 1, and 3 (see 
|26t [5^ I29j). Shown are the stability traces tr M(a) of two orbits P and F which are created from 
a period-tripling bifurcation of a straight-line orbit at a = 0.6315 and exist only for a < 0.6315. 
They are isolated at all values of a except for a = 0. At a = and near a ~ 0.0055 they appear to 
cross like in a TCB; note that also the same exchange of Maslov indices by one unit takes place for 
each orbit. 




-0.006 -0.003 0.0 0.003 0.006 0.009 

a 



Figure 12: Stability exchanges of orbits F and P in the coupled quartic oscillator. Both crossings 
are not transcritical bifurcations. (See text for details.) 

However, the two orbits have different shapes at all values of a, as shown by the inserts which 
display their shapes in the (x, y) plane. Therefore, their fixed points cannot coincide at either of 
the crossings, and hence the crossing near a ~ 0.0055 cannot be a TCB. The situation around 
a = 0, where the system is integrable, is not a bifurcation at all, but the generic Poincare-Birkhoff 
breaking of a rational (3:2) torus into a pair of stable and unstable isolated orbits. 

What actually happens near a ~ 0.0055, as described in [39j, is shown in Fig. [I3l where the 
scale of the graphs trM(a) has been enlarged in both directions by a factor ~ 10^. The graphs 
trM(a) really cross slightly above the critical line trM = 2 (note the shift along the vertical axis 
by two units), here at a distance of only 10^'^ which requires a high precision of the numerical 
calculations. Their crossings of the critical line, where bifurcations must take place, occurs at two 
different points in a tiny distance Aa = 3.56 x 10"''. These bifurcations are isochronous PFBs, and 
the orbit Q emitted and reabsorbed by them intermediates between the shapes of the two crossing 
orbits. Note that the orbit Q, which transforms the shape of the F orbit into that of the P orbit 
(or vice versa), has twice the discrete degeneracy of the orbits F and P. Compared to F this is 
because Q is a rotation and has two time orientations (while F transforms into itself under time 
reversal); compared to P it has a lower discrete symmetry (it is not symmetric under reflection at 
the X axis, but P is). This double degeneracy with respect to the parent orbits is characteristic of 
the isochronous (non-generic) PFB (see Sec. 12.2.2]) . 
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The change of the shape of the intermediate orbit Q is very similar to that of an orbit on an 
integrable torus - and can hardly be distinguished numerically from it. (It is exactly the case in 
the nearby integrable point a = 0, where P and F cross again and correspond to two realizations 
of one and the same 3:2 torus.) Under poor numerical resolution - as in Fig. [12] above - therefore, 
one might also misinterpret the situation at a ~ 0.0055 as existence of a locally integrable torus. 




a 



Figure 13: Same as Fig. [T2l but the scales in both directions are enlarged by a factor ~ 10^. Note 
how the orbit Q intermediates between the orbits F and P through two pitchfork bifurcations. 



Although this mechanism has been observed and published quite a long ago [39], it does not 
appear to be widely known. We deem it worth mentioning and illustrating here, in order to pre- 
vent misinterpretations of such "false" TCBs that can appear under poor numerical circumstances. 
(Although misunderstanding should not happen when the shapes of the orbits are known and/or 
when it is realized that their fixed points on the Poincare surface of section do not coincide at the 
crossing.) 
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